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Abstract 

The interaction between gravitational and electromagnetic radiation has a rather long 
research history. It is well known, in particular, that gravity-wave distortions can drive 
propagating electromagnetic signals. Since forced oscillations provide the natural stage for 
resonances to occur, gravito-electromagnetic resonances have been investigated as a means 
of more efficient gravity-wave detection methods. In this report, we consider the coupling 
between the Weyl and the Maxwell fields on a Minkowski background, which also applies 
to astrophysical environments where gravity is weak, at the second perturbative level. We 
use covariant methods that describe gravitational waves via the transverse component of the 
shear, instead of pure-tensor metric perturbations. The aim is to calculate the properties 
of the electromagnetic signal, which emerges from the interaction of its linear counterpart 
with an incoming gravitational wave. Our analysis shows how the wavelength and the am- 
plitude of the gravitationally driven electromagnetic wave vary with the initial conditions. 
More specifically, for certain initial data, the amplitude of the induced electromagnetic sig- 
nal is found to diverge. Analogous, diverging, gravito-electromagnetic resonances were also 
reported in cosmology. Given that, we extend our Minkowski-space study to cosmology 
and discuss analogies and differences in the physics and in the phenomenology of the Weyl- 
Maxwell coupling between the aforementioned two physical environments. 

PACS numbers: 98.80. Jk, 04.30.-w, 98.62.En 

1 Introduction 

Gravitational and electromagnetic waves are believed to interact in a variety of ways. The liter- 
ature contains a number of studies showing how gravitational radiation affects the propagation 
of electromagnetic signals by modifying their amplitude, their wavelength, their polarisation and 
their direction, either in vacuum or in the presence of conductive plasmas (e.g. see ^). There 
are also solutions of Einstein's equations that describe colliding plane gravitational and electro- 
magnetic waves on a flat Minkowski background. Only recently, however, techniques that could 
construct interacting solutions from the wave parameters before the collision were introduced 
(see ^ and references therein). Historically, most of the interest in the coupling between the 
Weyl and the Maxwell fields has been motivated by the on going search for gravitational-wave 
signals and by the quest for new and potentially more efficient detection mechanisms [3]. Several 
of the aforementioned references raise the possibility of resonances between the two sources and 
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then exploit it, either as a novel gravity- wave detection method [3], or as a general relativistic 
mechanism of amplifying large-scale magnetic fields 

Gravitational waves can affect electromagnetic signals directly (e.g. by modifying their am- 
plitude, wavelength, etc) as well as indirectly (e.g. by altering the geometry, or the nature, of the 
host environment) [S]. Here, we are interested in the direct effects, with particular emphasis on 
the possibility of resonances. The gravito-electromagnetic interaction is studied at second order 
on an empty, static Minkowski background. This means that our results apply to environments 
where the gravitational field is weak, like the interstellar space or the laboratory for example. In 
contrast to the majority of the approaches, which are metric based, we employ the 1-1-3 covariant 
approach to General Relativity |6] . The main difference is in the description of the gravitational 
radiation, which covariantly is monitored through the electric and magnetic components of the 
Weyl tensor [7]. The latter determines the free (i.e. the long range) gravitational field, the 
transverse part of which corresponds to gravity waves. 

Weyl curvature distortions affect the propagation of electromagnetic signals via Maxwell's 
equations. Gravitational waves trigger (pure-tensor) shear perturbations that directly affect 
both the electric and the magnetic component of the signal, through (the generalised) Ampere's 
and Faraday's laws. On a Minkowski background, the two sources propagate as simple (linear) 
plane waves, while their interaction is monitored by a second-order system that describes the 
gravitationally induced electromagnetic signal. Not surprisingly, the latter differs from the orig- 
inal one and the differences depend on the specifics of the gravito-electromagnetic interaction. 
For our purposes, however, the most important feature of the aforementioned second-order sys- 
tem is that it describes two forced oscillations (one for each component of the electromagnetic 
signal). The importance lies in the long known fact that forced vibrations provide the natural 
physical environment for resonances to occur. This means that, in principle at least, the inter- 
action between gravitational and electromagnetic waves can lead to the resonant amplification 
of the latter. 

Although resonances can occur spontaneously, they require rather special conditions. More 
specifically, the driving and the driven wave must oscillate in tune. In our case, the driving 
force of the gravitationally induced electromagnetic signal is described as a superposition of two 
waves. The characteristics of the latter, their effective wavelengths in particular, are fixed at the 
beginning of the interaction. The same initial conditions also determine the wavelength of the 
gravitationally induced - the driven - electromagnetic signal. We provide a set of formulae that 
allow one to calculate the specifics of the induced electromagnetic wave for a variety of initial 
conditions. We find that the wavelengths of the driving and the driven waves are generally dif- 
ferent, in which case the resulting electromagnetic oscillation will experience mild amplification 
(if any). When the circumstances are favourable, however, these two wavelengths will coincide 
and the driving and the driven source will oscillate in resonance. It is then, that the amplitude 
of the gravitationally induced electromagnetic signal is found to diverge. Whether this happens 
or not, according to our analysis, depends also on the interaction angle between the original elec- 
tromagnetic and gravitational waves. Here, we find that Weyl-Maxwell resonances occur when 
the aforementioned two interacting waves propagate in the same, or in the opposite, direction. 

These results should not come as a complete surprise. After all, we were dealing with forced 
oscillations, where resonances naturally occur. In fact, gravito-electromagnetic resonances have 
a rather long history in the relevant literature and certain types of them have been proposed in 
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the past as efficient mechanisms of gravitational-wave detection. To the best of our knowledge, 
however, the first time such strong (diverging) resonances were reported in the literature was in 
the context of cosmology Given that, and also for reasons of completeness and comparison, 
we extend our analysis to include cosmological environments as well. In cosmology, the host 
spacetime is no longer static, empty or flat. This means that, although the basic physics of 
the gravito-electromagnetic interaction remain essentially the same, the specifics change. The 
main difference appears to come from the expansion, which adds a "friction" term to the wave- 
equations that generally makes the resonance peak narrower. 

2 Electromagnetic and gravitational waves in Minkowski space 

In weak-gravity environments, for example in the interstellar space (away from massive compact 
stars) or in the laboratory, the spacetime is effectively flat. There, one can use the Minkowski 
metric to describe the geometry of the host environment. 

2.1 Electromagnetic waves 

Consider electromagnetic radiation traveling on an empty, static Minkowski background. The 
Maxwell field is invariantly described by means of the antisymmetric Faraday tensor (Fab), 
which splits into its electric and magnetic components once a family of timelike observers is 
introduced. Assuming that Ua is the 4-velocity of the aforementioned observers (normalised so 
that Uau"" = —1), the corresponding electric and magnetic fields are Ea = FabU^ and Ba = 
£abcF^'^/2 respectively {Sabc is the Levi-Civita tensor of the observers' 3-dimensional space). On 
our Minkowski background, the two fields obey the linear wave equations ^ 

Ea-B^Ea = and Ba-B^Ba = 0, (1) 

with the overdots indicating proper-time derivatives along the observers' worldlines. Also, = 
D'^Dtt represents the covariant Laplacian operating on the observers' 3-dimensional rest-space0 
Equations ([T^) and ([lb) are straightforward to solve. To proceed, we decompose the com- 
ponents of the electromagnetic signal by introducing the familiar harmonic splitting 

Ea = Y. Ein) and Ba = Y^ Qi2)(") , (2) 

n n 

where Qi^^^"\ qI^^^"^ are the corresponding vector harmonics and n is the eigenvalue/ wavenumber 
of the associated electromagnetic mode. Note that n? = n°'na, with representing the wavevec- 
tor along the propagation direction of the signal. Also, by construction, Da-E(„) = = Dai?(„), 

= = Qi'^^"\ D'^Qi^)^") = = D'^Qi')(") and D^.QiJ^^") = = D^^Sif ^"^1 Finally, we 
point out that both of the harmonic functions satisfy the vector version of the Laplace-Beltrami 
equation, that is 

D2g(l)(n) ^ _^2g(l)(n) D^Qi^)^ = -n^Qi^)^ . (3) 

^Given the flatness of the background metric, we have Da — >■ da and D'^ — D°Da — >■ d^da ~ to zero order. 
^Angled brackets indicate the symmetric and trace-free component of second-rank spacehke tensors. 
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Substituting decompositions ([2^) and ([2)3) back into Eqs. ([T^) and ([Da) and then using expres- 
sions dS]), the original wave formulae recast into 



and 



5(n) + n'B 



(n) 



0, 



(4) 



which now monitor the linear evolution of the signals' n-th harmonic mode. Both of the above 
are straightforward to solve, leading to the solution 



(n) 



£"0"^ sin (nto) H — Eq^^ cos (nto) 



+ 



Erf^ COS into) E""^^' sin (nto) 

n 



1 



sin {nt) 
cos (nt) , 



(5) 



for the electric component of the Maxwell field, and to solution 



B 



in] 



sin (nto) H — B^^' cos (nto) 



1 



n 



sin (nt) 



+ 



i^Q™'' COS (nto) — — Bq""' sin (nto) 



in), 



n 



cos (nt) 



(6) 



for its magnetic counterpart. This set monitors the propagation of an electromagnetic wave on 
an empty, static Minkowski background to first order. Note that the zero suffix indicates a given 
initial time, at which both the amplitude and the phase of signal are determined. 



2.2 Gravitational waves 

On the same Minkowski space, we will also consider the propagation of gravitational radiation. 
Covariantly, gravitational waves are described by means of the transverse parts of the electric 
and magnetic components of the Weyl field (e.g. see [7J). Nevertheless, given the symmetries of 
our background and the absence of matter, linear Weyl curvature distortions are monitored by 
the pure-tensor modes of the shear H These are isolated by imposing the first-order constraint 
D''(Jaft = (e.g. see [6]), which in our case is trivially satisfied at all times. Moreover, on a 
Minkowski background, the transverse component of the shear obeys the linear wave equation 

aab - DV„6 = . (7) 



As before, the overdots indicate time-differentiation relative to the fundamental observers and 
j-)2 _ represents their associated 3-D Laplacian operator. We also introduce the harmonic 

splitting 

k 

where (T(fc) is the A;- mode of the shear perturbation and Dacr(^k) = 0- Note that k'^ = k^ka, with k 
being the mode's eigenvalue/ wavenumber and ka the corresponding wavevector. Also, Q^^^ are 



(k) 



tensor harmonic functions, which satisfy the conditions Q^^ 
the tensor form 



Q 



(ab)' 



^ab 



= D^'Q^^^ and 



(k) 

ab 



-k'Q 



(fc) 



(9) 



'Also true for gravity waves propagating on a Friedmann- Robertson- Walker (FRW) background (see ^ 14. ip . 
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of the Laplace-Beltrami equation. Combining expressions ([7]), ([8]) and ([9|) leads to the hnear 
wave equation 



+ k (T(fc) 

which monitors the k-th harmonic mode of a gravitational wave propagating on an empty, static 
Minkowski background. As it is straightforward to show, the above accepts the solution 



0. 



(10) 



(k) 



a\^' sin (fcto) + t ct^q^ cos [kto] 



+ 



cos {kto) 



i a^Q^ sin (A;to) 



sin {kt) 
cos {kt) , 



(11) 



with the zero suffix corresponding to a given initial time and the quantities in brackets deter- 
mining the characteristics of the wave (i.e. its amplitude and phase). 



3 Gravito-electromagnetic interaction in Minkowski space 

Technically speaking, gravitational waves affect the propagation of their electromagnetic coun- 
terparts through Maxwell's equations. These combine to provide a set of two wave-equations 
(one for each component of the electromagnetic field), which allow us to study the coupling 
between the two sources at different perturbative levels. 



3.1 Second-order equations 

The last two sections considered the linear evolution of electromagnetic and gravitational ra- 
diation on an empty, static Minkowski background. Here, we will look into the nonlinear 
(second-order in particular) interaction of these two sources. Our aim is to study the effect 
of gravitational waves on electromagnetic signals. At the aforementioned perturbative level, the 
interaction between the Weyl and the Maxwell fields is governed by the system [Sj 

K - B^Ea = aabE + EabcBd^^'a'"' - 2eabc(r'dB^'B^^ - UabE' - EabE' + Hab& (12) 

and 

Ba - D^Ba = aabB - E abcE dl^^ cj"^ + 2e,b,CT^D<^^'^> - TZabB^ - Eab& - Hab& , (13) 

where Ea, Ba and Ba are the electromagnetic components before and after the interaction 
respectively. The former obey the simple, linear wave solutions ([5]) and ^ given in § 12.11 
Similarly, the transverse part of the shear satisfies solution (jlip . We also note that Rab is the 
perturbed (linear) 3-Ricci tensor, while Eab and Hab represent the (linearised on our Minkowski 
background) electric and magnetic parts of the Weyl tensor respectively. Then, to first order, 
we have [6] 

Tlab = Eab , Eab = -<^ab and Hab = cuilaab (14) 
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with curlua?) = £cd{a^'^'^'^b) by definition. Substituting the auxihary hnear relations (jl4p into 
Eqs. (fT2|) and p^ . the latter recast into 

- d2£;„ = aabE + 2aabE^ + eaftc^dOV^'^ - 2ea6et^^D<^^^> + B^curlcTab (15) 

and 

Ba - B^Ba = UabB + lUab^ " ^afec^dDV^'^ + 2e abc^^T)^' E'^'^ " ^^CUrl(T,b , (16) 

respectively. This set monitors the effects of gravitational radiation on electromagnetic signals 
propagating on a Minkowski background at second order. 

According to expressions (jl5p and (jl6p . the interaction between electromagnetic and gravi- 
tational radiation leads to electromagnetic waves that are generally different from their original 
counterparts. In particular, the specifics of the induced electromagnetic signals (i.e. their ampli- 
tude, wavelength and direction of propagation) depend on those of the interacting sources. For 
our purposes, however, the most important implication of (jlSp . ()16p is that gravitational- wave 
distortions can drive electromagnetic signals. This is important, because forced oscillations pro- 
vide the natural physical environment for resonances to occur. In other words, Eqs. (jl5p and 
(|16p imply that the interaction between the Weyl and the Maxwell field can lead to the resonant 
amplification of the latter. 

3.2 Resonant couplings 

Resonances require rather special conditions, which makes them more likely to occur in controlled 
environments than spontaneously in nature. If gravito-electromagnetic resonances are possible, 
the system of (fT5]l . (fTB]) should contain resonant solutions. To demonstrate (analytically) that 
this is indeed the case, let us ignore the backreaction of the electric component upon its magnetic 
counterpart and vice versa. Then, Eqs. (fT5]) and (fT6]l reduce to 

Ea - T^^Ea = aabE + 2&ab& (17) 

and ^ 

Ba-D^Ba = aabB +2&abB\ (18) 

respectively. To simplify the mathematics, without compromising the basic physics, suppose 
that both the electromagnetic signal and the gravitational wave are monochromatic. Then, we 
may use the decompositions 

^a = %)Qi')(") and ^, = %)Qi2)W, (19) 

for the components of the original electromagnetic field and 

aab = a^k)Qil\ (20) 

for the gravitationally induced shear perturbation. As mentioned before, the characteristics of 
the gravitationally driven (the induced) electromagnetic wave depend on those of the originally 
interacting sources. This dependance is encoded in the decomposition 

Ea = E^e) Qi'^^') and B^ = qI^) W , (21) 
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where Qa = Sit Q(i)(ra) Qa = Q(2)(n) by construction. These are vector harmonic 
functions, with qI^^^^^ = = Qi^^^^\ D-^qI^^^^^ = = B^Q^a^^^l Also, 

D^qWW = and D^Qi^)^ ^ _^2g(2)(£) ^ (22) 

with 

£2 = „2 ^ /c2 _^ 2nk cos . (23) 

The last expression relates the wavelength of the gravitationally driven electromagnetic wave to 
the wavelengths and the interaction angle (</>, with < < vr) of the original sources. More 
specifically, given the values of n, k and </>, one can employ (j23p to obtain a unique value for i. 

Substituting decompositions ([T2D - (^1]) back into expression (fT7|) and using the Laplace- 
Beltrami equation (|22fe i). we arrive at 

% + e^E^e) = <y{k)E + 2a(fc)^(") . (24) 

The driving force in the right-hand side is determined by the interaction between the original 
electromagnetic signal and the incoming gravitational wave. Following the discussion given in 
§ 12.11 and § 12. 2( these are monitored by 

= ^sin(nt + i?^) and (T(fc) = Bsui{kt + i?e) , (25) 

with the amplitudes and the phases depending on the initial conditions (see solutions ([5]) and 
(fTTI) ). Without loss of generality, we may set the phases in (f25k ) and (|25b ) to zero. Then, 
expression (pi]) leads to 

+ fE(^i,) = ^C{n- 2k) sin [{k -n)t] + ^ C{n + 2k) sin [{k + n)t] , (26) 

where C = AB. Clearly, proceeding in an exactly analogous way, we obtain the wave equation 
of the gravitationally driven magnetic component 

+ fB(^(,-) = ^C{n- 2k) sin [{k - n) t] + ^ C{n + 2k) sin [{k + n)t\ . (27) 

The last two equations show that the force driving the propagation of the gravitationally induced 
electromagnetic signal can be expressed as the superposition of two waves, the specifics of which 
are determined by the initial conditions. In particular, the effective wavenumbers of the driving 
oscillations are mi^2 = A:=pn, while that of the gravitationally driven electromagnetic signal {() is 
given by (j23p . Clearly, when the initial conditions are favourable, these wavenumbers coincide. 
Then, both the driving and the driven waves will oscillate in resonance. 

3.3 Resonant solutions 

The simplest resonant case occurs for k = n, that is when the original sources share the same 
wavelength (though their directions of propagation may differ). Then, it is straightforward to 
show that the pair of (p6]) . (p7|) reduces to 

+ fEi^g) = Tsm{mt) and + fB^e^ = J'sin(mt) , (28) 
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with T = 'iCk/2 and m = 2k. In other words, the driving force can be expressed as a single 
wave with amphtude and wavelength depending on those of the original sources. The solutions 
of ([^5h ) and ([25b ) are respectively given by 

E^i-j = W sin {it + (/?) + -g2^-^ sin (mt) (29) 

and 

J" 

= W sin (£t + (^) + ^2_^2 ^™ ("^*) ' (20) 

where the values of the integration constants (W and if) are fixed at the onset of the gravito- 
electromagnetic interaction. Both of the above show the amplitude of the gravitationally driven 
(the induced) electromagnetic wave to diverge as £ — > m. According to the relation (j23p . this 
corresponds to — )• 0. Consequently, the interaction between electromagnetic and gravitational 
radiation, which share the same wavelength and propagate in the same direction on a Minkowski 
background, will resonantly amplify the electromagnetic signal. 

Generally, the wavelengths of the interacting gravitational and electromagnetic signals will 
be different (i.e. k ^ n). More specifically, in most physically realistic situations we expect 
that the wavelength of the gravitational radiation will far exceed that of its electromagnetic 
counterpart. This implies that typically n ^ k. In addition, the two originally interacting 
waves will generally have nonzero (and in principle different) phases. When k ^ n and t?^, 
i?e 7^ 0, a similar analysis shows that the solution of Eq. (j24p takes the form 



7"^ J^2 

= ^ sin {It + 93) + 2 + '^i) + 72 2 ^'^^ ("^2* + (^2) , (31) 

where J-"i^2 = C{n^2k) /2, mi^2 = k^n and ^1^2 = ''^A-f^B- It goes without saying that an exactly 
analogous solution describes the propagation of the gravitationally driven magnetic component 
as well. Following (j3ip . the amplitude of the electromagnetic wave diverges as £ — > |?7ii,2|- 
In the first case this translates to i — ?■ |fc — n|, which occurs when (p ^ tt (see Eq. (I23p ). 
When i k + n, on the other hand, we find that (p 0. Therefore, the interaction between 
electromagnetic and gravitational waves, which have different wavelengths but propagate in 
the same or in the opposite direction on a Minkowski background, will resonantly amplify the 
electromagnetic signal. 

4 Gravito-electromagnetic interaction in cosmology 

The basic physics leading to the Weyl-Maxwell resonances of the previous section are generic 
and not particular to the geometrical structure of the Minkowski space. This means that, at 
least in principle, analogous effects should also take place within more general spacetimes, like 
those related to cosmology for example. 

4.1 Background and linear dynamics 

In cosmology, the spacetime is no longer empty, static and fiat, which means that the specifics 
of the interaction between the Weyl and the Maxwell fields generally change. Here, we will 
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consider the gravito-electromagnetic coupling on an FRW background with Euchdean spatial 
hypersurfaces. We will assume, in particular, a spatially flat Friedmann universe containing a 
single barotropic fluid. The time evolution of this model is determined by the zero-order set 

H'^ = -p, H = -H^ -- p(l + 3w) and p = -3Hp(l + w), (32) 
3 6 

where H = a /a is the Hubble parameter (a is the cosmological scale factor), p is the density 
and w = p/p is the barotropic index of the matter (with p representing the isotropic pressure). 

Let us now perturb the above defined background with propagating gravitational and elec- 
tromagnetic waves0 As before, the former are monitored by the transverse component of the 
shear, which now obeys the linear wave equation (e.g. see [6]) 

(Tab - D^0"a6 = -5H&ab - -H'^{1- 3w)aab ■ (33) 

At the same perturbative level, and given that our FRW background has flat spatial sections, 
the magnetic field decays adiabatically at all times. This means that, to linear order, 

Ba oc , (34) 

irrespective of the type of matter that fills the universe and of its electric properties (i.e. of 
whether the electrical conductivity is high or low [9j). Next, we will look into the interaction 
between these two sources at the second perturbative level. 

4.2 Second order gravito-magnetic interaction 

Inflation is probably the best available mechanism capable of causally producing large-scale 
gravitational and electromagnetic perturbations. These distortions start as subhorizon quantum 
fluctuations and freeze-out as classical perturbations once outside the Hubble scale (e.g. see [lOj). 
Also, throughout the de Sitter phase, the cosmic medium is believed to be a very poor electrical 
conductor. After inflation, during reheating and the subsequent radiation era, the electrical 
conductivity of the matter increases rapidly. As a result, the electric fields gradually vanish and 
the currents freeze the magnetic field in with the fluid. Put another way, the post-inflationary 
universe satisfies the ideal-magnetohydronamic (MHD) requirements. Nevertheless, causality 
ensures that the MHD limit holds only well within the Hubble scale. Beyond the horizon, 
however, the electric currents had no time to establish themselves and the low-conductivity 
approximation of the de Sitter period still holds. 

In what follows, we will consider the gravito-magnetic interaction at the time the Weyl and 
the Maxwell fields start crossing back inside the Hubble radius. Close to the horizon, the electric 
currents are not yet sufficiently strong to freeze the magnetic fields in with the matter. Thus, the 
electric fields are still present, though weakened by the ever increasing strength of the currents. 

''The gravito-magnetic interaction was considered in [4], assuming a weakly magnetised FRW background. 
This meant that, strictly speaking, the mechanism was linear in perturbative terms. Here, to facilitate a direct 
comparison with the Minkowski-space study of the previous section, we adopt a magnetic-free background and 
consider the coupling between the Weyl and the Maxwell fields at the second perturbative level. 
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All these mean that, near the Hubble threshold, the effect of the gravitational waves on the 
magnetic fields is monitored by the second-order expression [8] 

B^Ba + 5HBa + 3(1 - w)H^Ba = 2 i&ah + 2Haab) & • 



Ba 



(35) 



In analogy with the Minkowski-space case discussed earlier, Ba represents the gravitationally 
induced magnetic field. The main difference between the above and Eq. (|18|) in ^ 13.21 is in the 
expansion and the matter terms seen in the left-hand side of ([35]). As in the Minkowski case, the 
most important feature (for our purposes) of the latter expression is the gravito-magnetic term 
in its right-hand side. The presence of this particular source-term means that Weyl-curvature 
distortions can drive magnetic oscillations in cosmological environments as well. In other words, 
Weyl-Maxwell resonances can also occur between cosmological gravitational waves and large- 
scale, primordial magnetic fields. 



4.3 Cosmological gravito-magnetic resonances 

The outcome of the gravito-magnetic interaction is determined by the coupled system of ()33p - 
(|35]) and by the evolution of our background model (see Eqs. (p2]l ). To proceed further, we 
harmonically decompose the perturbations in a manner analogous to that used in § 13.21 Then, 
expressions ([33l) and (i35]l take the form 



2 ^ 



2,w)H' 



'k 
.a 



a, 



(k) 







and 



B{e) + 5HB(i) + 



(^) 



3{1 -w)H^ + 



B. 



(36) 



(37) 



, we have = B^^^{ao/aY 



respectively. Note that i is the eigenvalue/ wavenumber of the gravitationally driven i?-field, k is 
that of the gravitational wave, n determines the scale of the background i?-field and the three are 
related by Eq. (j23p [3]. Also, according to the linear relation ([3l 
at all times, where the zero suffix specifies a moment during the universe's evolution. 

Current observations indicate that the current Hubble length lies close to 10^ Mpc. This 
implies that most of the astrophysically relevant scales have entered the horizon during the late 
radiation era. Then, solving (p6|) for w = 1/3, a oc t^^'^ and H = l/2t, we obtain 



(fc) 



k 



cos 



k 



aoHo 



1 




(38) 



with 6 = tan~^{aQ&Q^^ /kuo'^) representing the phase and the zero suffix indicating the beginning 
of the radiation epoch. Substituting the above into the right-hand side of Eq. (I37p . while taking 
into account that -B(„) oc always, we arrive at 



2t2 



+ 



-)■ 

00/ 



B, 



2B, 



(n) 



k\a); 



ao 



5/2 



X sm 



aoHo 




(39) 
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The latter monitors the gravitationally induce magnetic field as it crosses inside the Hubble 
horizon of a radiation dominated, spatially flat FRW universe. Like its Minkowski counterpart 
before, expression (j39p describes a forced oscillation, where the driving force on the right-hand 
side is the result of the gravito-magnetic interaction. Unlike its Miknowski analogue, on the 
other hand, Eq. (|39|) contains extra terms due to the expansion of the universe and the presence 
of matter. Another difference is in the amplitude of the driving wave, which is no longer constant 
but decays in time. 

The nature of p9p means that resonant solutions, analogous to those found in Minkowski- 
space, are in principle possible in cosmology as well. The additional (expansion-related) terms, 
on the other hand, suggest that the specifics of the resonances may differ. Indeed, solving 
Eq. ([Ml) we obtain (see also § IV.A and Eq. (22) in [4]) 



B 



B, 



in) 



±- 



X sm 



iHo \tj MaoFoVio, 




^i"Vfc2(a(^))^ + ag(a(^))^ fto\ (iTk [T\ 
^Ho [tj \aoHo\toJ 



X cos 




(40) 



where Si(x) and Ci{x) represent the sine and the cosine integral functions respectively [llj. Fol- 
lowing the above, the gravitationally induced S-field generally oscillates with an amplitude that 
decays adiabatically in time (i.e. -B(f) oc oc a~^). At any given time, however, the amplitude 
depends on the initial conditions. This dependence is crucial because it also determines the ar- 
gument of the aforementioned integral functions. The latter are known to satisfy the conditions 
lim3;^oSi(x) = and limi;_>.o Ci(x) = — oo [11]. Therefore, the presence of the cosine integral 
function in solution (140 p . ensures that the amplitude of the gravitationally driven magnetic field 
diverges when the two sources oscillate in resonance (i.e. for k ^ tj. The same integral functions 
also mean that resonances between the Weyl and the Maxwell fields in cosmology have a much 
narrower peak than their Minkowski space analogues. This is due to the universal expansion, 
which acts as a dumping force and weakens the overall effect of the gravito-electromagnetic 
resonance. Broader resonance peaks can be achieved, when the source term in the right-hand 
side of Eq. ([39]) drops slower with time. In particular, the resonant amplification of the induced 
-B-field is exactly analogous to that seen in Eq. (I30p . when the amplitude of the driving wave 
depletes as instead of This can happen, for example, when the background magnetic 

field decays as (i.e. B oc t~^/'^ during the radiation era), like many of the i?-fields proposed 
in various scenarios of primordial magnetogenesis [9j. Then, one obtains 



2aoi?^"Vfc2(a, 
fc2 



B 



^0 . 



sm 



k 

aoHf) 




(41) 
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for the gravitationally driven magnetic field [1]. The above solution shows resonant magnetic 
amphfication, when k ^ £, analogous to that of the Minkowksi background (compare (j4ip to 
Eq. dSni) earlier). 

Note that, in principle at least, the amplification effect of the resonances seen in solutions 
(j4U|) and (j41|) is independent of the amount of energy stored in gravity-wave perturbations. 
In practise, this means that even the weak gravitational waves that inflation is expected to 
produced, can resonantly amplify primordial (most likely also inflationary) magnetic fields during 
the subsequent evolution of the universe. We should also point out that the amplification 
mechanism discussed above operates around the time of the second horizon crossing. There the 
accompanying primordial electric fields are still present, though weakened by the ever growing 
effect of the currents. Deep inside the Hubble radius, however, the currents will dominate, 
eliminate the electric fields and freeze their magnetic counterparts into the highly conductive 
cosmic fluid. Once the freezing-in process is completed, the ideal MHD limit will be established. 
Then onwards, resonances are no longer possible and our amplification mechanism ceases to 
operate [4J. In the mean time, however, the universe has been permeated by large-scale magnetic 
fields that are substantially stronger than the original ones and, perhaps, strong enough to seed 
the galactic dynamo [12j . 

5 Discussion 

The basic physical principles behind the gravito-electromagnetic interaction are simple and 
essentially independent of the environment. More specifically, gravitational-wave distortions in- 
teract and can act as driving forces of electromagnetic radiation, irrespective of the structure of 
the host spacetime. This means that gravito-electromagnetic resonances are in principle possible 
both in astrophysics and cosmology, as well as in the laboratory. Nevertheless, there are still 
major differences between these environments, which are expected to affect the phenomenology 
of the resulting resonances. Among the goals of the present study was to identify such analo- 
gies and differences. The approach we have adopted describes gravitational-wave perturbations 
via the transverse component of the shear, rather than by using the pure-tensor part of the 
metric perturbation. This facilitates a direct and physically more transparent coupling between 
the gravitational and the electromagnetic fields through Maxwell's equations. In line with the 
purposes of our study, the Weyl-Maxwell coupling was considered first on an empty Minkowski 
background and the results were then compared to those obtained within a perturbed FRW 
cosmology. In either case, we have allowed for the presence of gravitational and electromagnetic 
waves and analysed the effects of the former on the latter at second order. In both cases, we 
have found that the system of equations monitoring the gravito-electromagnetic interaction con- 
tained resonant solutions. To be precise, under favourable initial conditions - when the driving 
and the driven waves oscillate in tune, the amplitude of the gravitationally induced electromag- 
netic/magnetic field was found to diverge. There were, however, differences as well. Applied to 
Minkowski space, the gravito-electromagnetic interaction led to a textbook-like wave equation 
that describes forced oscillations free of any dumping terms. In cosmology, however, the expan- 
sion of the universe acts as a source of dumping that weakens the resonance, in the same (more 
or less) way that inhibits the growth of density perturbations. This weakening is refiected in 
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the resonance peak, which in cosmology is generally narrower. 

Resonances can occur either spontaneously in nature, or under controlled conditions in the 
laboratory. In either case, weak vibrations lead to oscillations of disproportionately large ampli- 
tude. It is not surprising therefore that resonances between gravitational and electromagnetic 
waves have been discussed as a potential gravity-wave detection mechanism. This way, the 
ubiquitous presence of electromagnetic radiation could prove a very useful ally in our ongoing 
effort to detect the still elusive gravity- wave signals. Cosmology is also an area where resonances 
between the Weyl and the Maxwell fields could have an effect, by amplifying weak primordial 
magnetic fields to astrophysically relevant strengths. Then, one may be able to address the still 
open issue of cosmic magnetogenesis, within the realm of standard cosmology and conventional 
electrodynamics. In principle, analogous gravito-electromagnetic resonances should also occur 
on spacetimes with more complicated structures, though then one would probably have to use 
numerical (rather than analytical) methods to identify them. Overall, in view of the generic 
nature of the Weyl-Maxwell coupling within general relativity and given that forced oscillations 
provide the natural stage for resonances to occur, the question may be how often such resonances 
take place, rather than whether they happen at all. This remains to be seen. If the answer turns 
out to be "positive" , however, the ubiquitous presence of magnetic and electromagnetic fields in 
our cosmos may has been partly achieved at the expense of gravitational radiation. 
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Abstract 

The interaction between gravitational and electromagnetic radiation has a rather long 
research history. It is well known, in particular, that gravity-wave distortions can drive 
propagating electromagnetic signals. Since forced oscillations provide the natural stage for 
resonances to occur, gravito-electromagnetic resonances have been investigated as a means 
of more efficient gravity-wave detection methods. In this report, we consider the coupling 
between the Weyl and the Maxwell fields on a Minkowski background, which also applies 
to astrophysical environments where gravity is weak, at the second perturbative level. We 
use covariant methods that describe gravitational waves via the transverse component of the 
shear, instead of pure-tensor metric perturbations. The aim is to calculate the properties of 
the electromagnetic signal, which emerges from the interaction of its linear counterpart with 
an incoming gravitational wave. Our analysis shows how the wavelength and the amplitude 
of the gravitationally driven electromagnetic wave vary with the initial conditions. More 
specifically, for certain initial data, the amplitude of the induced electromagnetic signal 
is found to diverge. Analogous, diverging, gravito-electromagnetic resonances were also 
reported in cosmology. Given that, we extend our Minkowski-space study to cosmology 
and discuss analogies and differences in the physics and in the phenomenology of the Weyl- 
Maxwell coupling between the aforementioned two physical environments. 

1 Introduction 

Gravitational and electromagnetic waves are believed to interact in a variety of ways. The liter- 
ature contains a number of studies showing how gravitational radiation affects the propagation 
of electromagnetic signals by modifying their amplitude, their wavelength, their polarisation and 
their direction, either in vacuum or in the presence of conductive plasmas (e.g. see [T]). There 
are also solutions of Einstein's equations that describe colliding plane gravitational and electro- 
magnetic waves on a flat Minkowski background. Only recently, however, techniques that could 
construct interacting solutions from the wave parameters before the collision were introduced 
(see [2J and references therein). Historically, most of the interest in the coupling between the 
Weyl and the Maxwell fields has been motivated by the on going search for gravitational-wave 
signals and by the quest for new and potentially more efficient detection mechanisms [3]. Several 
of the aforementioned references raise the possibility of resonances between the two sources and 
then exploit it, either as a novel gravity- wave detection method [3j, or as a general relativistic 
mechanism of amplifying large-scale magnetic fields [1]. 



1 



Gravitational waves can affect electromagnetic signals directly (e.g. by modifying their am- 
plitude, wavelength, etc) as well as indirectly (e.g. by altering the geometry, or the nature, of the 
host environment) [5]. Here, we are interested in the direct effects, with particular emphasis on 
the possibility of resonances. The gravito-electromagnetic interaction is studied at second order 
on an empty, static Minkowski background. This means that our results apply to environments 
where the gravitational field is weak, like the interstellar space or the laboratory for example. In 
contrast to the majority of the approaches, which are metric based, we employ the 1+3 covariant 
approach to General Relativity [6j . The main difference is in the description of the gravitational 
radiation, which covariantly is monitored through the electric and magnetic components of the 
Weyl tensor [7J. The latter determines the free (i.e. the long range) gravitational field, the 
transverse part of which corresponds to gravity waves. 

Weyl curvature distortions affect the propagation of electromagnetic signals via Maxwell's 
equations. Gravitational waves trigger (pure-tensor) shear perturbations that directly affect 
both the electric and the magnetic component of the signal, through (the generalised) Ampere's 
and Faraday's laws. On a Minkowski background, the two sources propagate as simple (linear) 
plane waves, while their interaction is monitored by a second-order system that describes the 
gravitationally induced electromagnetic signal. Not surprisingly, the latter differs from the orig- 
inal one and the differences depend on the specifics of the gravito-electromagnetic interaction. 
For our purposes, however, the most important feature of the aforementioned second-order sys- 
tem is that it describes two forced oscillations (one for each component of the electromagnetic 
signal). The importance lies in the long known fact that forced vibrations provide the natural 
physical environment for resonances to occur. This means that, in principle at least, the inter- 
action between gravitational and electromagnetic waves can lead to the resonant amplification 
of the latter. 

Although resonances can occur spontaneously, they require rather special conditions. More 
specifically, the driving and the driven wave must oscillate in tune. In our case, the driving 
force of the gravitationally induced electromagnetic signal is described as a superposition of two 
waves. The characteristics of the latter, their effective wavelengths in particular, are fixed at the 
beginning of the interaction. The same initial conditions also determine the wavelength of the 
gravitationally induced - the driven - electromagnetic signal. We provide a set of formulae that 
allow one to calculate the specifics of the induced electromagnetic wave for a variety of initial 
conditions. We find that the wavelengths of the driving and the driven waves are generally dif- 
ferent, in which case the resulting electromagnetic oscillation will experience mild amplification 
(if any). When the circumstances are favourable, however, these two wavelengths will coincide 
and the driving and the driven source will oscillate in resonance. It is then, that the amplitude 
of the gravitationally induced electromagnetic signal is found to diverge. Whether this happens 
or not, according to our analysis, depends also on the interaction angle between the original elec- 
tromagnetic and gravitational waves. Here, we find that Weyl-Maxwell resonances occur when 
the aforementioned two interacting waves propagate in the same, or in the opposite, direction. 

These results should not come as a complete surprise. After all, we were dealing with forced 
oscillations, where resonances naturally occur. In fact, gravito-electromagnetic resonances have 
a rather long history in the relevant literature and certain types of them have been proposed in 
the past as efficient mechanisms of gravitational-wave detection. To the best of our knowledge, 
however, the first time such strong (diverging) resonances were reported in the literature was in 
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the context of cosmology [3]. Given that, and also for reasons of completeness and comparison, 
we extend our analysis to include cosmological environments as well. In cosmology, the host 
spacetime is no longer static, empty or flat. This means that, although the basic physics of 
the gravito-electromagnetic interaction remain essentially the same, the specifics change. The 
main difference appears to come from the expansion, which adds a "friction" term to the wave- 
equations that generally makes the resonance peak narrower. 

2 Electromagnetic and gravitational waves in Minkowski space 

In weak-gravity environments, for example in the interstellar space (away from massive compact 
stars) or in the laboratory, the spacetime is effectively flat. There, one can use the Minkowski 
metric to describe the geometry of the host environment. 

2.1 Electromagnetic waves 

Consider electromagnetic radiation traveling on an empty, static Minkowski background. The 
Maxwell field is invariantly described by means of the antisymmetric Faraday tensor (Fab), 
which splits into its electric and magnetic components once a family of timelike observers is 
introduced. Assuming that Ua is the 4-velocity of the aforementioned observers (normalised so 
that Mqu" = —1), the corresponding electric and magnetic fields are Ea = Fahv!' and Ba = 
EabcF^'^/'^ respectively {Sabc is the Levi-Civita tensor of the observers' 3-dimensional space). On 
our Minkowski background, the two fields obey the linear wave equations [8] 

Ea-B'^Ea = and Ba-B^Ba = 0, (1) 

with the overdots indicating proper-time derivatives along the observers' worldlines. Also, = 
D^Dfl represents the covariant Laplacian operating on the observers' 3-dimensional rest- spaceS 
Equations ([I^) and ([U)) are straightforward to solve. To proceed, we decompose the com- 
ponents of the electromagnetic signal by introducing the familiar harmonic splitting 

K = E ^(") and Ba = Y: B^n) Qi')^") , (2) 

n n 

where Qi^^^"^ qI^^^"^ are the corresponding vector harmonics and n is the eigenvalue/ wavenumber 
of the associated electromagnetic mode. Note that ii? = n°'na, with representing the wavevec- 
tor along the propagation direction of the signal. Also, by construction, Dai?(„) = = Dai?(„), 
q(i)(") = = Q(2)(n)^ D«Q«(") = = D"Qi')(") and D^.Q^J^^") = = D^.sg^^"). Finally, we 
point out that both of the harmonic functions satisfy the vector version of the Laplace-Beltrami 
equation, that is 

D2g{l)(n) ^ _^2g(l)(n) D^Qi^)^ = -n^Qi^)^ . (3) 

Substituting decompositions ([2^) and dSJa) back into Eqs. ([1^) and ([U)) and then using expression 
([3]), the original wave formulae recast into 

+ n2^(„) = and + n^B^^^ = , (4) 

Given the flatness of the background metric, we have T)a -> da and = D"Da d^da = 9^ to zero order. 
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which now monitor the hnear evolution of the signals' n-th harmonic mode. Both of the above 
are straightforward to solve, leading to the solution 



(n) 



Ejp^ sin (nto) H — -Eq"^ ("-^o) 



n 



sin (nt) 



E^^ cos (nto) - — Eq^' sin (nto) 



n 



cos (nt) , 



(5) 



for the electric component of the Maxwell field, and to solution 



B 



(n) 



B^^ sin (nto) H — B^""' cos (nto) 



1 



>(") 



n 



sin (nt) 



+ 



i^Q™'* COS (nto) Bq""' sin (nto) 



in), 



n 



cos (nt) 



(6) 



for its magnetic counterpart. This set monitors the propagation of an electromagnetic wave on 
an empty, static Minkowski background to first order. Note that the zero suffix indicates a given 
initial time, at which both the amplitude and the phase of signal are determined. 



2.2 Gravitational waves 

On the same Minkowski space, we will also consider the propagation of gravitational radiation. 
Covariantly, gravitational waves are described by means of the transverse parts of the electric 
and magnetic components of the Weyl field (e.g. see [7]). Nevertheless, given the symmetries of 
our background and the absence of matter, linear Weyl curvature distortions are monitored by 
the pure-tensor modes of the shearH These are isolated by imposing the first-order constraint 
D"(Ta5 = (e.g. see [6j), which in our case is trivially satisfied at all times. Moreover, on a 
Minkowski background, the transverse component of the shear obeys the linear wave equation 



0. 



(7) 



As before, the overdots indicate time-differentiation relative to the aforementioned observers and 
j-)2 _ represent their associated 3-D Laplacian operator. We also introduce the harmonic 

splitting 

^ab = E^W2i?' (8) 

k 

where cr(fc) is the fe-mode of the shear perturbation and Da(T(fc) = 0. Note that k"^ = k^ka, with 
k being the mode's eigenvalue/ wavenumber and ka the corresponding wavevector. Also, Q^^J 



(k) 



are tensor harmonic functions, which satisfy the conditions Q^^^ 
the tensor form 



2(S>> Safe = = D''q1J and 

(9) 



of the Laplace-Beltrami equation. Combining expressions ([7]), ([8]) and ([9|) leads to the linear 
wave equation 



(k) 



0. 



(10) 



This is also true for gravitational waves propagating on an Friedmann-Robertson- Walker (FRW) background. 
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which monitors the A;-th harmonic mode of a gravitational wave propagating on an empty, static 
Minkowski background. As it is straightforward to show, the above accepts the solution 



0"i 



(fe) 



alP^ sin (kto) + — dg"''' cos {kto] 
K 



An) 



+ 



Uq""^ cos {ktQ) 



(Tq""^ sin (kto) 



k 



sin [kt) 
cos {kt) , 



(11) 



with the zero sufhx corresponding to a given initial time and the quantities in brackets deter- 
mining the characteristics of the wave (i.e. its amplitude and phase). 



3 Gravito-electromagnetic interaction in Minkowski space 

Theoretically, gravitational waves affect the propagation of the electromagnetic counterparts 
through Maxwell's equations. These combine to provide a set of two wave-equations (one for 
each component of the electromagnetic field), which allow us to study the coupling between the 
two sources at different perturbative levels. 



3.1 Second-order equations 

The last two sections considered the linear evolution of electromagnetic and gravitational radia- 
tion on an empty, static Minkowski space. Here, we will look into the nonlinear (second-order in 
particular) interaction of these two sources. Our aim is to study the effect of gravitational waves 
on electromagnetic signals. At the aforementioned perturbative level, the interaction between 
the Weyl and the Maxwell fields is governed by the system [8] 

Ea - = aahE + e^bc^dDV^'^ - 2e,b,a^D<^5'^) - Uab^ - EabE^ + HabB^ (12) 

and 

Ba - B^Ba = aabB - eabcEdT^'(7''' + 2eabc(T'dB^'E^^ - UabB' - Eab& - HabE' , (13) 

where Ea, Ba and Ea, Ba are the electromagnetic components before and after the interaction 
respectively. The former obey the simple, linear wave solutions ([5]) and ([6j) given in § 12.11 
Similarly, the transverse part of the shear satisfies solution (jlip . We also note that Rat is the 
perturbed (linear) 3-Ricci tensor, while Eab and Hab represent the (linearised on our Minkowski 
background) electric and magnetic parts of the Weyl tensor respectively. Then, to first order, 
we have [6] 

T^ab = Eab , Eab = -(^ab and Hab = CUrliTab (14) 

with curler^;, = ecd(aD'^^'^b) by definition. Substituting the auxiliary linear relations p4p into 
Eqs. (fT2]) and (fT3|) . the latter recast into 

Ea - B^Ea = aabE + 2aab& + eabcBaB^'a''^ - 2e„b,a^D<^^'^) + ^^curla„6 (15) 
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and 

Ba - B'^Ba = cJabB + 2dab& " e^fec^dDV^'^ + 2e,b,cT^D<^^'^> - &cm\aab , (16) 

respectively. This set monitors the effects of gravitational radiation on electromagnetic signals 
propagating on a Minkowski background at second order. 

According to expressions (jl5|) and (|16|) . the interaction between electromagnetic and gravi- 
tational radiation leads to electromagnetic waves that are generally different from their original 
counterparts. In particular, the specifics of the induced electromagnetic signals (i.e. their ampli- 
tude, wavelength and direction of propagation) depend on those of the interacting sources. For 
our purposes, however, the most important implication of (llSp . (116p is that gravitational- wave 
distortions can drive electromagnetic signals. This is important, because forced oscillations pro- 
vide the natural physical environment for resonances to occur. In other words, Eqs. (jlSp and 
(jl6p imply that the interaction between the Weyl and the Maxwell field can lead to the resonant 
amplification of the latter. 

3.2 Resonant couplings 

Resonances require rather special conditions, which makes them more likely to occur in controlled 
environments than spontaneously in nature. If gravito-electromagnetic resonances are possible, 
the system of (llSp . (I16p should contain resonant solutions. To demonstrate (analytically) that 
this is indeed the case, let us ignore the backreaction of the electric component upon its magnetic 
counterpart and vice versa. Then, Eqs. (fT5]) and (fT6]l reduce to 

Ea - Ti^Ea = UabE + Id^bE^ (17) 

and ^ 

Ba-B^Ba = aabB +2&abB\ (18) 

respectively. To simplify the mathematics, without compromising the physics, suppose that 
both the electromagnetic signal and the gravitational wave are monochromatic. Then, we may 
use the decompositions 

^a = %)Qi')(") and ^, = %)Qi2)(n)^ (19) 

for the components of the original electromagnetic field and 

^afe = ^(fc)Qi?, (20) 

for the gravitationally induced shear perturbation. As mentioned before, the characteristics of 
the gravitationally driven (the induced) electromagnetic wave depend on those of the originally 
interacting sources. This dependance is encoded in the decomposition 

Ea = E^e)Q^^^^'^ and = i3(,)Q(2)W , (21) 

where q2'^^^^ = Q^ab ^\i){n) ^^'^ Q^a^^^^ = 2ife'*2(2)(n) by Construction. These are vector harmonic 
functions, with qI^^^^^ = = Qi^^^^\ D'^qI^^^"^ = = D^'qI^^^^^ Also, 

D^qWW = -£2q{i)W and D^Q^W = ^ (22) 
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with 

f =n'^ -^k'^ + 2nk cos (j) . (23) 

The last expression relates the wavelength of the gravitationally driven electromagnetic wave to 
the wavelengths and the interaction angle ((/>, with < < vr) of the original sources. More 
specifically, given the values of n, k and (p, one can employ (|23|) to obtain a unique value for i. 

Substituting decompositions ([T9|) - ([2T|) back into expression (fTT]) and using the Laplace- 
Beltrami equation (122b ,). we arrive at 

% + fE^,) = + 2a(fc)^W . (24) 

The driving force in the right-hand side is determined by the interaction between the original 
electromagnetic signal and the incoming gravitational wave. Following the discussion given in 
§ 12.11 and § 12. 2| these are monitored by 

= ^sin(nt + and (T(fc) = Bsm{kt + i?^) , (25) 

with the amplitudes and the phases depending on the initial conditions (see solutions ([5]) and 
PT|) ). Without loss of generality, we may set the phases in ([^5h ) and (|25b ) to zero. Then, 
expression (pi]) leads to 

1 1 

+ = -C{n- 2k) sin [{k -n)t] + - C{n + 2k) sin [{k + n)t] , (26) 

where C = AB. Clearly, proceeding in an exactly analogous way, we obtain the wave equation 
of the gravitationally driven magnetic component 

B(^i) + fB(^i) = i C(n - 2k) sin [{k - n) t] + ^ C{n + 2k) sin [{k + n)t] . (27) 

The last two equations show that the force driving the propagation of the gravitationally induced 
electromagnetic signal can be expressed as the superposition of two waves, the specifics of which 
are determined by the initial conditions. In particular, the effective wavenumbers of the driving 
oscillations are mi^2 = k^n, while that of the gravitationally driven electromagnetic signal (i) is 
given by l\23h . Clearly, when the initial conditions are favourable, these wavenumbers coincide. 
Then, both the driving and the driven waves will oscillate in resonance. 

3.3 Resonant solutions 

The simplest resonant case occurs for k = n, that is when the original sources share the same 
wavelength (though their directions of propagation may differ). Then, it is straightforward to 
show that the pair of (f26]) . (i27|l reduces to 

+ = J'cos(mt) and + fB(^i) = J'cos(mt) , (28) 

with = 3Ck/2 and m = 2k. In other words, the driving force can be expressed as a single 
wave with amplitude and wavelength depending on those of the original sources. The solutions 
of (f28k ) and (i28b) are respectively given by 

E(^e) = W sin {et + ip)+ 2 5 2 ("^*) (29) 
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and 

Ba) = W sin {it + ip)+ „^ . sin {mt) , (30) 

where the value of integration constants (W and ip) are fixed at the onset of the gravito- 
electromagnetic interaction. Both of the above show the amphtude of the gravitationally driven 
(the induced) electromagnetic wave to diverge as £ m. According to the relation (j23b ). this 
corresponds to — ?> 0. Consequently, the interaction between electromagnetic and gravitational 
radiation, which share the same wavelength and propagate in the same direction on a Minkowski 
background, will resonantly amplify the electromagnetic signal. 

Generally, the wavelengths of the interacting gravitational and electromagnetic signals will 
be different (i.e. k ^ n). More specifically, in most physically realistic situations we expect 
that the wavelength of the gravitational radiation will far exceed that of its electromagnetic 
counterpart. This implies that typically n ^ k. In addition, the two originally interacting 
waves will generally have nonzero (and in principle different) phases. When k ^ n and t?^, 
T?e 7^ 0, a similar analysis shows that the solution of Eq. (j24p takes the form 

= W sin {£t + if) + 2 ^™ ("^1* + ^1) + 2 ("^2* + 0J2) , (31) 

where J-"i^2 = C(n=p2fc)/2, mi^2 = k^n and uji^2 = ^A-f'^B- It goes without saying that an exactly 
analogous solution describes the propagation of the gravitationally driven magnetic component 
as well. Following (j3ip . the amplitude of the electromagnetic wave diverges as ^ — > |mi^2|- 
In the first case this translates to £ — )• \k — n\, which occurs when (/>—)• tt (see Eq. (I23|) ). 
When £ — > fc + n, on the other hand, we find that (/> — > 0. Therefore, the interaction between 
electromagnetic and gravitational waves, which have different wavelengths but propagate in 
the same or in the opposite direction on a Minkowski background, will resonantly amplify the 
electromagnetic signal. 



4 Gravito- electromagnetic interaction in cosmology 

The basic physics leading to the Weyl-Maxwell resonances of the previous section are generic 
and not particular to the geometrical structure of the Minkowski space. This means that, at 
least in principle, analogous effects should also take place within more general spacetimes, like 
those related to cosmology for example. 



4.1 Background and linear dynamics 

In cosmology, the spacetime is no longer empty, static and flat, which means that the specifics of 
the interaction between the Weyl and the Maxwell fields generally change. Here, we will consider 
the gravito-electromagnetic coupling on an Friedmann-Robertson- Walker (FRW) background 
with Euclidean spatial hyper surf aces. We will assume, in particular, a spatially flat Friedmann 
universe containing a single barotropic fluid. The time evolution of this model is determined by 
the zero-order set 

if2 = -p, H = -H^ -- p(l + 3w) and p = -3Hpn+w), (32) 
3 6 
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where H = a /a is the Hubble parameter (a is the cosmological scale factor), p is the density 
and w = p/p is the barotropic index of the matter (with p representing the isotropic pressure). 

Let us now perturb the above defined background with propagating gravitational and elec- 
tromagnetic wavesH As before, the former are monitored by the transverse component of the 
shear, which now obeys the linear wave equation (e.g. see [B]) 

3 

(Jab - D^CJah = -bH&ab - - - 3w)(Taf, . (33) 

At the same perturbative level, and given that our FRW background has flat spatial sections, 
the magnetic field decays adiabatically at all times. This means that, to linear order, 

Ba oc , (34) 

irrespective of the type of matter that fills the universe and of its electric properties (i.e. of 
whether the electrical conductivity is high or low [9]). Next, we will look into the interaction 
between these two sources at the second perturbative level. 

4.2 Second order gravito-magnetic interaction 

Inflation is probably the best available mechanism capable of causally producing large-scale 
gravitational and electromagnetic perturbations. These distortions start as subhorizon quantum 
fluctuations and freeze-out as classical perturbations once outside the Hubble scale (e.g. see pi)]). 
Also, throughout the de Sitter phase, the cosmic medium is believed to be a very poor electrical 
conductor. After inflation, during reheating and the subsequent radiation era, the electrical 
conductivity of the matter increases rapidly. As a result, the electric fields gradually vanish and 
the currents freeze the magnetic field in with the fiuid. Put another way, the post-inflationary 
universe satisfies the ideal-magnetohydronamic (MHD) requirements. Nevertheless, causality 
ensures that the MHD limit holds only well within the Hubble scale. Beyond the horizon, 
however, the electric currents had no time to establish themselves and the low-conductivity 
approximation of the de Sitter period still holds. 

In what follows, we will consider the gravito-magnetic interaction at the time the Weyl and 
the Maxwell fields start crossing back inside the Hubble radius. Close to the horizon, the electric 
currents are not yet sufficiently strong to freeze the magnetic fields in with the matter. Thus, the 
electric fields are still present, though weakened by the ever increasing strength of the currents. 
All these mean that, near the Hubble threshold, the effect of the gravitational waves on the 
magnetic fields is monitored by the second-order expression [8] 

Ba - B^Ba + 5HBa + 3(1 - w)H^Ba = 2 i&ab + 2Haab) B^ ■ (35) 

In analogy with the Minkowski-space case discussed earlier, Ba represents the gravitationally 
induced magnetic field. The main difference between the above and Eq. (llSp in § 13.21 is in the 

^The gravito-magnetic interaction was considered in [4|, assuming a weakly magnetised FRW background. 
This meant that, strictly speaking, the mechanism was linear in perturbative terms. Here, to facilitate a direct 
comparison with the Minkowski-space study of the previous section, we adopt a magnetic-free background and 
consider the coupling between the Weyl and the Maxwell fields at the second perturbative level. 
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expansion and the matter terms seen in the left-hand side of (j35p . As in the Minkowski case, the 
most important feature (for our purposes) of the latter expression is the gravito-magnetic term 
in its right-hand side. The presence of this particular source-term means that Weyl-curvature 
distortions can drive magnetic oscillations in cosmological environments as well. In other words, 
Weyl-Maxwell resonances can also occur between cosmological gravitational waves and large- 
scale, primordial magnetic fields. 



4.3 Cosmological gravito-magnetic resonances 



The outcome of the gravito-magnetic interaction is determined by the coupled system of 
(j35]) and by the evolution of our background model (see Eqs. ([32]) ). To proceed further, we 
harmonically decompose the perturbations in a manner analogous to that used in § 13. 2[ Then, 
expressions ([33|) and ([35]) take the form 



1 (1 - 3w)H' + 

2 \a 



cr. 



(k) 







and 



3(1 -w)H^ + ' ^ 



B. 



(36) 



(37) 



respectively. Note that i is the eigenvalue/ wavenumber of the gravitationally driven 5-field, k is 
that of the gravitational wave, n determines the scale of the background 5-field and the three are 
related by Eq. (j23|) [4j. Also, according to the linear relation (p^ . we have = B^^^{ao/a)'^ 
at all times, where the zero suffix specifies a moment during the universe's evolution. 

Current observations indicate that the current Hubble length lies close to 10^ Mpc. This 
implies that most of the astrophysically relevant scales have entered the horizon during the late 
radiation era. Then, solving (|36p for w = 1/3, a cx and H = l/2t, we obtain 



(fc) 



k [j ' 




(38) 



with 6 = taW^laoUQ^^ /kuQ^^) representing the phase and the zero suffix indicating the beginning 
of the radiation epoch. Substituting the above into the right-hand side of Eq. (i37|) . while taking 
into account that B^^-j oc always, we arrive at 



2t2 + UoJ U 



B, 




X sm 



5/2 



(39) 



The latter monitors the gravitationally induce magnetic field as it crosses inside the Hubble 
horizon of a radiation dominated, spatially flat FRW universe. Like its Minkowski counterpart 
before, expression (j39p describes a forced oscillation, where the driving force on the right-hand 
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side is the result of the gravito-magnetic interaction. Unhke its Miknowski analogue, on the 
other hand, Eq. (j39p contains extra terms due to the expansion of the universe and the presence 
of matter. Another difference is in the amplitude of the driving wave, which is no longer constant 
but decays in time. 

The nature of (|39p means that resonant solutions, analogous to those found in Minkowski- 
space, are in principle possible in cosmology as well. The additional (expansion-related) terms, 
on the other hand, suggest that the specifics of the resonances may differ. Indeed, solving 
Eq. dMI) we obtain (see also § IV.A and Eq. (22) in g]) 



B 



{(■) 



±- 




X sm 




X cos 




(40) 



where Si(x) and Ci(a;) represent the sine and the cosine integral functions respectively pTj. Fol- 
lowing the above, the gravitationally induced S-field generally oscillates with an amplitude that 
decays adiabatically in time (i.e. oc cx a~^). At any given time, however, the amplitude 
depends on the initial conditions. This dependence is crucial because it also determines the ar- 
gument of the aforementioned integral functions. The latter are known to satisfy the conditions 
lim2,._5.o Si(x) = and lim^_^o Ci(x) = — oo [H]. Therefore, the presence of the cosine integral 
function in solution ()40p . ensures that the amplitude of the gravitationally driven magnetic field 
diverges when the two sources oscillate in resonance (i.e. for k i). The same integral functions 
also mean that resonances between the Weyl and the Maxwell fields in cosmology have a much 
narrower peak than their Minkowski space analogues. This is due to the universal expansion, 
which acts as a dumping source and weakens the overall effect of the gravito-electromagnetic 
resonance. Broader resonance peaks can be achieved, when the source term in the right-hand 
side of Eq. (I39p drops slower with time. In particular, the resonant amplification of the induced 
B-Held is exactly analogous to that seen in Eq. (jSOp . when the amplitude of the driving wave 
depletes as instead of This can happen, for example, when the background magnetic 

field decays as (i.e. B gc during the radiation era), like many of the i?-fields proposed 
in various scenarios of primordial magnetogenesis [9]. Then, one obtains 



2aoi?S"V^'(^i 



B 



^0 . 



— sm 




(41) 



for the gravitationally driven magnetic field [3]. The above solution shows resonant magnetic 
amplification, when k ^ analogous to that of the Minkowksi background (compare (j4ip to 
Eq. ([30]) earlier). 
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Note that, in principle at least, the amplification effect of the resonances seen in solutions 
(j40p and (j4ip is independent of the amount of energy stored in gravity-wave perturbations. 
In practise, this means that even the weak gravitational waves that inflation is expected to 
produced, can resonantly amplify primordial (most likely also inflationary) magnetic fields during 
the subsequent evolution of the universe. We should also point out that the amplification 
mechanism discussed above operates around the time of the second horizon crossing. There the 
accompanying primordial electric fields are still present, though weakened by the ever growing 
effect of the currents. Deep inside the Hubble radius, however, the currents will dominate, 
eliminate the electric fields and freeze their magnetic counterparts into the highly conductive 
cosmic fluid. Once the freezing-in process is completed, the ideal MHD limit will be established. 
Then onwards, resonances are no longer possible and our amplification mechanism ceases to 
operate [1]. In the mean time, however, the universe has been permeated by large-scale magnetic 
fields that are substantially stronger than the original ones and, perhaps, strong enough to seed 
the galactic dynamo [12] . 

5 Discussion 

The basic physical principles behind the gravito-electromagnetic interaction are simple and 
essentially independent of the environment. More specifically, gravitational-wave distortions in- 
teract and can act as driving forces of electromagnetic radiation, irrespective of the structure of 
the host spacetime. This means that gravito-electromagnetic resonances are in principle possible 
both in astrophysics and cosmology, as well as in the laboratory. Nevertheless, there are still 
major differences between these environments, which are expected to affect the phenomenology 
of the resulting resonances. Among the goals of the present study was to identify such analogies 
and differences. The approach we have adopted in this study, describes gravitational-wave per- 
turbations via the transverse component of the shear, rather than by using the pure-tensor part 
of the metric perturbation. This facilitates a direct and physically more transparent coupling 
between the gravitational and the electromagnetic fields through Maxwell's equations. In line 
with the purposes of our study, the Weyl-Maxwell coupling was considered first on an empty 
Minkowski background and the results were then compared to those obtained within a perturbed 
FRW cosmology. In either case, we have allowed for the presence of gravitational and electro- 
magnetic waves and analysed the effects of the former on the latter at second order. In both 
cases, we have found that the system of equations monitoring the gravito-electromagnetic inter- 
action contained resonant solutions. To be precise, under favourable initial conditions - when 
the driving and the driven waves oscillate in tune, the amplitude of the gravitationally induced 
electromagnetic/magnetic field was found to diverge. There were, however, differences as well. 
Applied to Minkowski space, the gravito-electromagnetic interaction led to a textbook-like wave 
equation that describes forced oscillations free of any dumping terms. In cosmology, however, 
the expansion of the universe acts as a source of dumping that weakens the resonance, in the 
same (more or less) way that inhibits the growth of density perturbations. This weakening is 
reflected in the resonance peak, which in cosmology is generally narrower. 

Resonances can occur either spontaneously in nature, or under controlled conditions in the 
laboratory. In either case, weak vibrations lead to oscillations of disproportionately large ampli- 
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tude. It is not surprising therefore that resonances between gravitational and electromagnetic 
waves have been discussed as a potential gravity-wave detection mechanism. This way, the 
ubiquitous presence of electromagnetic radiation could prove a very useful ally in our ongoing 
effort to detect the still elusive gravity- wave signals. Cosmology is also an area where resonances 
between the Weyl and the Maxwell fields could have an effect, by amplifying weak primordial 
magnetic fields to astrophysically relevant strengths. Then, one may be able to address the still 
open issue of cosmic magnetogenesis, within the realm of standard cosmology and conventional 
electrodynamics. In principle, analogous gravito-electromagnetic resonances should also occur 
on spacetimes with more complicated structures, though then one would probably have to use 
numerical (rather than analytical) methods to identify them. Overall, in view of the generic 
nature of the Weyl-Maxwell coupling within general relativity and given that forced oscillations 
provide the natural stage for resonances to occur, the question may be how often such resonances 
take place, rather than whether they happen at all. This remains to be seen. If the answer turns 
out to be "positive" , however, the ubiquitous presence of magnetic and electromagnetic fields in 
our cosmos may has been partly achieved at the expense of gravitational radiation. 
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